function [ ] = ctrPlotsWrapper(ctrM, ctrC, ets, edat, figures_path, optfig, step_func_macs, originalPriceModel)

%% Plot all step function MAC curves versus iso-elastic
%if nargin > 6 && step_func_macs == true
%    filestub  = [ figures_path '/fig_step_mac_pd' ];
%    plotAllStepFuncMACs(ctrM, filestub, optfig );
%    close all;
%end

%% Aggregate MAC Curve
%filestub = [ figures_path '/fig_agg_mac' ];
%plotAggCostFuncWrapper(ctrM, filestub, optfig)

% Time series of prices by compliance period
filename = [ figures_path '/fig_time_series_price.pdf' ];
if nargin > 7
    plotFitTimeSeriesPrice( ctrM.output_tk, edat.period, filename, optfig, originalPriceModel.output_tk );
else
    plotFitTimeSeriesPrice( ctrM.output_tk, edat.period, filename, optfig);
end


%% Distribution of emissions, data and model
filestub  = [ figures_path '/fig_hist_emissions_' ];
titlestub = 'Distribution of emissions';
inlogs    = true;
plotFitHistEmissionsAll( edat.plant, ctrM, edat.period, ...
    filestub, titlestub, inlogs, optfig );
close all;


%% Scatter emissions in simulation against data (model generated) by period
filestub  = [ figures_path '/fig_scatter_emissions' ];
titlestub = 'Scatter plot of emissions in model vs. data';
inlogs    = true;
plotFitScatterEmissionsAll( edat.plant, ctrM, edat.period, ...
    filestub, titlestub, inlogs, optfig )
close all;

% Should this comparison be using a different emissions series?
%   Bidding data series may include validations not present in the
%   emissions data as constructed.


%% Run plot comparing costs in market and command-and-control regimes
titlestub = 'Variable abatement cost by regime';
filestub = [ figures_path '/fig_total_cost_comparison_' ];
plotTotalCostsAll( ctrM, ctrC, edat.period,...
    filestub, titlestub, optfig )
close all;

% TK: Check on treatment of non-operational plants
end
